We present a new density-matrix functional within the 
recently introduced framework for tensor-product expansions 
of the two-particle density matrix. It performs well both for 
the homogeneous electron gas as well as atoms. For the ho- 
mogeneous electron gas, it performs significantly better than 
all previous density-matrix functionals, becoming very accu- 

jm -«. . rate for high densities and outperforming Hartree-Fock at 

metallic valence electron densities. For isolated atoms and 

ions, it is on a par with generalized gradient approximations 

to density-functional theory. We also present analytic re- 

*T~^' suits for the correlation energy in the low density limit of 

r j. . , the free electron gas for a broad class of such functionals. 
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The long standing quest for an accurate description of 
electronic correlation with effective single-particle theories 
has renewed interest in density-matrix functionals in recent 
years. The success of density functional approaches, based on 
Hohenberg-Kohn-Sham theory MM, is well established |3J. 
Unfortunately, the known approximations to the requisite uni- 
versal energy functional of the density are not sufficiently ac- 
curate to predict the rates of chemical and physical processes 
at room temperatures. This theory can easily be generalized 
to show that there is a universal correlation energy functional 
of the single-particle density matrix as well H-ftNIM . A den- 
sity matrix functional requires an expression of the two-body 
matrix in terms of the one-body density matrix. An ansatz 
for the two-body matrix containing one free parameter was 
proposed by Miiller JlQ| , and the case where this parameter 
is set to 1/2 has been discussed by Buijse fill]. In addition, 
Buijse attempted to perform self-consistent calculations but 
optimized only the lowest two natural orbitals and occupa- 
tion numbers and only for the H2 molecule. Independently, 
Goedecker and Umrigar |l2J proposed a functional based on a 
very similar, but self-interaction corrected, two-body matrix 
and derived the Euler-Lagrange formalism necessary for the 
numerical minimization of the functional. Optimizing typi- 
cally 50 natural orbitals and occupation numbers, necessary 
to achieve convergence, for a variety of atoms and ions it was 
found that the Goedecker-Umrigar (GU) functional yielded 
energies that were comparable or better than those from the 
generalized gradient approximation (GGA) in density func- 
tional theory and densities that were better than those from 
GGA. Following these first encouraging numerical results sev- 
eral groups of authors have recently proposed and studied 
the accuracy of various new approximations to the density- 
matrix functional |13-b(|. However, numerical studies fl4f] 
demonstrate that density-matrix functionals, to date, give 
poor correlation energies for the homogeneous electron gas. 

It has been shown analytically that the per particle corre- 
lation energy of the homogeneous electron ga s in the low den- 
sity limit approaches a non-zero constant [ff5j for these initial 
functionals. Building on our tensor-product formalism [ J14J , 
below we extend this proof and demonstrate this shortcoming 
for a broad class of functionals. As regions of low densities 
contribute relatively little to the total energies of solids and 
molecules, we focus our attention on improving the results 
for the high and intermediate density regimes. We present a 
new density-matrix functional which is more than an order of 
magnitude more accurate in the high-density limit for the ho- 
mogeneous electron gas than those proposed previously and 
also represents a significant improvement in the intermediate 
regime. In atoms, the new functional performs similarly to 
the generalized gradient approximation (GGA) fell to den- 
sity functional theory but not as well as the GU functional. 



Because selfinteraction corrections complicate numerical cal- 
culations, we do not enforce them on our new functional in 
the present work. 

Notation — For a system of N electrons, the exact energy 
functional of the iuio-body density matrix 7(xiX2,x' 1 x 2 ) = 



E = / dxi 



-i\7 2 ri +U( ri ) 



n(xi,xi; 



(1) 



+ / dxidx; 



7(xiX 2 ,XiX 2 ) 

|ri -r 2 | 



Here and throughout, we work in atomic units, and the Xi 
are compound coordinates representing both position rj and 
spin Si, so that integration over X; represents integration over 
space and summation over spin channels. In the above ex- 
pression, the external potential and one-body density matrix 
are U(r) and n(xi,x' 1 ) = (^|^ t (xi)i/>(xi)|^'}, respectively. 
Finally, we note that simple integration of the definitions re- 
veals the sum-rule which connects the one- and two-body 
matrices, 



( — - — Jn(xi,xi)= / dx2 7(xiX2,xix2). 



(2) 



As in (14|, we expand 7 into Hartree and exchange- 
correlation contributions, 



7 = 7h + 7xc- 



(3) 



By separation of variables, it is always possible to expand the 
unknown exchange-correlation part as a (potentially infinite) 
sum of tensor products, 



7xc = 2J «» ® & 



(4) 



In this context we view the two-body density matrix as a two- 
body operator (a function of four variables) and Ui and vt as 
one-body operators (functions of two variables). Note that 
throughout we employ the twisted tensor product (denoted 
as "Type III" product in Q), 

[u®u](xixi,X2X 2 ) = w(xi,X 2 )t>(xi,X 2 ), 

to combine the one-body operators. 

Remarkably, given such an expansion truncated to finite or- 
der, the combination of the symmetry constraints on 7 with 
the sum-rule is sometimes sufficient to determine the full 
four-variable function 7 directly in terms of the two-variable 
density matrix n. Under these conditions, one can then eval- 
uate the energy as an explicit functional of the density matrix 
using (Jlh . Minimizing the density-matrix functional over all 



physical density matrices then yields the ground-state of the 
system under this approximation. 

New Functional — Considerable freedom remains in the 
terms for the expansion (Mj). Previously, we explored two 
choices, "Corrected Hartree-Fock" (CHF) Q, and the "Cor- 
rected Hartree" (CH) functional, 



-CHF 

7xc 



-CH 

7xc 
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-(y/njl^h)' 



i/n(l - n)), 



respectively. The latter is identical to the one proposed 
by Buijse [hl| and a special case of the one proposed by 
Miiller (hJ. The GU functional is the CH functional with 
a selfinteraction correction. In the case of the homogeneous 
electron gas at high densities the GU and CH functionals be- 
come identical since the selfinteraction correction vanishes but 
at low densities they could differ if the natural orbitals are lo- 
calized. One way to view these functionals is in terms of the 
coefficient of the Fock term n <g> n, which is —1/2 in the CHF 
functional and is zero in the CH functional. Under truncation 
of the expansion to one additional term, the value of this co- 
efficient, through the sum-rule, then uniquely determines the 
form of the remaining tensor product. 




FIG. 1. Correlation energies of the three functionals in the 
text and the exact result J23-E4J, versus the mean particle 
separation. The solid lines are numerical calculations, the 
dashed lines are exact results. The vertical lines indicate the 
range of typical metallic valence densities. 

Homogeneous electron gas — Figure [D shows the correlation 
energies of the resulting functionals for the spin-unpolarized 
homogeneous electron gas. Corrected Hartree-Fock leads to 
no improvement in the correlation energy over Hartree-Fock, 
and approaches zero for high electron density. The CH func- 
tional is a little better, and represents an improvement by 
a factor of two over Hartree-Fock at high densities, where 
it is over-correlated. Both functionals are over-correlated at 
low densities, where they perform badly, with the per parti- 
cle correlation energy approaching a non-zero constant. This 
property of the CH functional has been demonstrated analyt- 
ically in Q. 

Having seen the behaviour of the above two functionals in 
the high density regime, it is natural to look for a form with 



a Fock coefficient intermediate to those above. A coefficient 
of —1/4 for the Fock term results in the following functional, 



7xc 



-- (h® h + y/h(2 - h) \/h(2- n)j 



(5) 



Figure hi shows that this new functional results in a dramatic 
improvement at high densities (low r s ) and significant im- 
provement at metallic valence densities (1 < r s < 6). At low 
density, the behaviour of all three functionals becomes the 
same, resulting in unrealistically large correlation energies per 
particle. Note that because this happens at low density, the 
contribution of these regions to the total energy is likely to be 
small in an inhomogeneous system. We now explore the rea- 
sons for this behaviour to ascertain the feasibility of making 
improvements in this regime as well. 

Low density limit — To understand the approach of the 
three functionals shown in Fig. fil to one another in the low 
density limit, we recall that Cioslowski and Pernal [J15| have 
determined the predictions of the CH functional analytically 
for r s > 5.769. Their analysis applies, provided that 7xc is 
represented as a single tensor product of a specific homoge- 
neous form. As an illustration of the anomalous correlation 
behaviour which these functionals exhibit in the low density 
limit, Figure compares the known many-body momentum 
distribution with numerical results for our new functional and 
the analytic distribution from ]l5[ derived for the CH func- 
tional for r s = 10. Both density-matrix functionals have very 
low occupation numbers, and lose the characteristic signature 
of the sharp drop near the Fermi wave vector evident in the 
exact curve. 
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FIG. 2. The momentum distribution }(k) as a function of 
the wave vector, k: exact |25j] (dashed line), CH functional 
[ p.5| (dash-dotted line) and our new functional (solid line) at 
r 3 = 10. The Fermi wave vector is denoted fc/ 

We now show that, as a consequence of the sum-rule, all 
of the above functionals reduce to this particular form in the 
limit of high r s . Fourier transforming (hi) , the total energy in 
the homogeneous electron gas becomes 

'V^dkdk' 



E 



Vdk k 2 



/k 



4tt 



-ffkk' 



(6) 



(2tt) 3 2 J V J (2tt) 6 ,,, ,,, 

where <? kk / is the momentum space representation of 7xc- In 
the case that 7xc is expanded as a sum of tensor products 



of single particle operators, its momentum space representa- 
tion will be a sum of tensor products of the momentum space 
representation of those operators. Because our one-body op- 
erators are explicit functions of the density matrix, which 
for the homogeneous electron gas is diagonal in the momen- 
tum representation, their momentum space representations 
are also diagonal. Thus, their tensor products become simple 
products, and we have 



3kk' 



fii[n(k), n(k')] + fi 2 [n(k), n(k')] + 



(7) 



where the fl are the same functions of the density ma- 
trix which appear in the tensor products expressing 7xc 
but with the tensor-product operator replaced with simple 
multiplication. For example, for the CH functional <? kk / = 
-y n (k) n (k'). 

We now show that, in the low-density limit, the analytic 
solution presented in |l5|| for the case when g kk i consists of a 
single Q, function of the homogeneous form 



fi[An(k), An(k)] = A /3 n[n(k),n(k')], 



(8) 



applies more generally to cases with multiple tensor products. 
The analysis of [[15| begins with the scaling ansatz 

n(k) = p A iq{p"\s.), 

where p is the total electron density per spin channel, v is a 
constant, r\ is a universal function that does not depend on 
the density explicitly, and A = 1 — 3v due to the constraint 
of the total particle density. Note that as p — > 0, n(k) — * 
uniformly (provided that A > 0). This behaviour, evident in 
Figure g, forms the crux of our argument; we verify it self- 
consistently below. 

In the above limit, all of the Qj in the above functionals ap- 
proach homogeneous functions with their own exponents j3i. 
Substituting n(k) into (m), we then find that a single homo- 
geneous function dominates g kk ', the one with the smallest 
Pi (= /3 m in). Therefore, all the conditions for the analysis of 
5[ apply, in particular that 
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One can show that not only the above functionals, but 
any functional satisfying the sum-rule and dominated by a 
single homogeneous term in the low-density limit, will exhibit 
Anin = 1 in that limit. If indeed a single tensor product term 
dominates 7, then the symmetry arguments presented in ]14] 
imply that 7xc must take the form 

7xc = u[n] ® fi[fi]. (9) 

The sum-rule (g) then implies 

rfx^xc = n, (10) 

which, using (B), can be written as an operator equation 

(u[h]) 2 =h, (11) 

which implies that u[ ] must be a homogeneous function with 
exponent 1/2 and therefore that 7xc is dominated by a term 



of homogeneous exponent /3 m m = 1. Therefore v — and 
A = 1; the fact that A > confirms the self-consistency and 
the vali dity of the above analysis. 

Ref. fllq] analyzes the case f3 = 1 in detail, finding for the 
correlation energy per particle the result 



as p — > 0. Figure nl verifies that, in the low density limit, 
indeed the correlation energy of all three functionals approach 
one another and, in particular, that they all approach a non- 
zero constant correlation energy of —1/8 Hartree per particle. 

Atoms — The electrons in solids exhibit both free-electron 
like and localized atomic-orbital behaviour. Thus, another 
important limiting case for assessing new functionals is their 
behavior in atoms. Table | presents energies obtained from 
the new functional for light atoms and ions. For comparison, 
we include energies obtained from the CH and GU functionals 
as well as from the local density approximation (LDA), and 
the Perdew-Burke-Ernzerhof (PBE) pi| generalized gradient 
approximation. 

The table reveals that, on the whole, the new functional 
performs for atoms quite comparably to both the CH func- 
tional and to PBE. In particular, for the two-electron se- 
ries, the new functional is somewhat worse, but for the four- 
electron series, it is somewhat better. All density-matrix 
functionals are generally much better than LDA. For the neu- 
tral atoms, CH and the present functional are not quite as 
good as PBE, but, for the ions both are better or in some 
cases similar to PBE. Note that the CH functional always 
overestimates the correlation energy, which is not the case for 
the new functional. In the GU functional this overcorrelation 
is corrected by the selfinteraction correction and this func- 
tional yields the best energies for atoms and ions. Thus, we 
feel that the density-matrix functionals provide an adequate 
description of atom-like systems. 

Conclusions — We have introduced a new density-matrix 
functional, which represents a major improvement over exist- 
ing functionals of the density matrix. It is very accurate in 
the high density regime of the homogeneous electron gas, sig- 
nificantly improves the correlation energies at typical valence 
densities, and it is comparable to the generalized gradient 
approximation in atoms. Although it is not as accurate as 
the GU functional for atoms, we expect it to perform bet- 
ter than the GU functional for solids. We have also shown 
that, in our tensor-product expansion of the two-body den- 
sity matrix, little further improvement can be expected at the 
present truncation of two terms. Our test systems span the 
range of environments encountered in solids and molecules, 
so we conclude that this new functional is a good candidate 
to be used in electronic structure calculations of condensed 
matter. 

We thank Cyrus Umrigar for interesting discussions and 
helpful comments on the manuscript. 
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Atom 


Present 


CH 


GU 


LDA 


PBE 


He 


0.01 


-0.01 


0.006 


0.07 


0.01 


Be 


-0.07 


-0.1 


-0.02 


0.2 


0.04 


Ne 


0.1 


-0.07 


0.05 


0.7 


0.07 


Be 2 + 


-0.02 


-0.005 


0.004 


0.2 


0.04 


C 4+ 


-0.02 


-0.003 


0.003 


0.4 


0.06 


C 2+ 


0.1 


-0.2 


0.01 


0.4 


0.07 


Q 4+ 


0.1 


-0.2 


0.02 


0.6 


0.1 



TABLE I. Error in the total energies of atoms for the new 
functional, in Hartrees. For comparison, results are included 
for the CH and the GU functionals and for the LDA and 
PBE-GGA approximations within density functional theory. 



